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Abstract 

This work aims at providing a mathematical and numerical framework for the analysis on the 
effects of pulsed electric fields on biological media. Biological tissues and cell suspensions are described 
as having a heterogeneous permittivity and a heterogeneous conductivity. Well-posedness of the model 
problem and the regularity of its solution are established. A fully discrete finite element scheme is 
proposed for the numerical approximation of the potential distribution as a function of time and space 
simultaneously for an arbitrary shaped pulse, and it is demonstrated to enjoy the optimal convergence 
order in both space and time. The presented results and numerical scheme have potential applications 
in the fields of medicine, food sciences, and biotechnology. 
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1 Introduction 

The electrical properties of biological tissues and cell suspensions determine the pathways of current 
flow through the medium and, thus, are very important in the analysis of a wide range of biomedical 
applications and in food sciences and biotechnology mmm- 

A biological tissue is described as having a permittivity and a conductivity m- The conductivity 
can be regarded as a measure of the ability of its charge to be transported throughout its volume by an 
applied electric field while the permittivity is a measure of the ability of its dipoles to rotate or its charge 
to be stored by an applied external field. At low frequencies, biological tissues behave like a conductor 
but capacitive effects become important at higher frequencies due to the membranous structures [22] [23]. 

Iu this paper, we consider a model problem for the effect of pulsed electric fields on biological tissues. 
Our goal is to study the electric behavior of a biological tissue under the influence of a pulsed electric 
field. It is of great importance to understand the effects of the pulse shape on the potential distribution 
in the tissue medium. We provide a numerical scheme for computing the potential distribution as a 
function of time and space simultaneously for an arbitrary shaped pulse. Our results are expected to 
have important applications in neural activation during deep brain simulations 0111 , debacterization of 
liquids, food processing [23], and biofouling prevention [2Tj. Our numerical scheme can be also used for 
selective spectroscopic imaging of the electrical properties of biological media [2]. It is challenging to 
specify the pulse shape in order to give rise to selective imaging of cell suspensions mm- 

The paper is organized as follows. In section [2] we introduce the model equation and some notations 
and preliminary results. We recall the method of continuity and the notions of weak and strong solutions. 
Section [3] is devoted to existence, uniqueness, and regularity results for the solution to the model problem. 
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We first derive an a priori energy estimate. Then we prove existence and uniqueness of the weak solution. 
Finally, we investigate the interface problem where the conductivity and permittivity distributions may 
be discontinuous, which is a common feature of biological media. It is shown in section [3] that the 
solution to the interface problem has a higher regularity in each individual region than in the entire 
domain. This regularity result is critical for our further numerical analysis. In section [I] we investigate 
the numerical approximation of the solution to the interface problem. Assuming that the domain is a 
convex polygon, we present a semi-discrete scheme and prove the error estimates for it in both H 1 - and 
L 2 -norms. With these estimates at hand, we then process to propose a fully-discrete scheme and establish 
the error estimates for it in both H 1 - and L 2 -norms. It is worth mentioning that both semi-discrete and 
fully discrete scheme achieve optimal convergence order in both H 1 - and L-norm, provided that the 
interface is exactly resolved. 

Let us end this section with some notations used in this paper. For a domain U C R ra ,n = 2 or 
3, each integer k > 0 and real p with 1 < p < oo, W k,p (U) denotes the standard Sobolev space of 
functions with their weak derivatives of order up to k in the Lebesgue space L P {U). When p = 2, we 
write H k (U) for W k,2 {U). The scalar product of L 2 (fl) is denoted by (•,•). If A is a Banach space 
with norm || • || A and J C K is an interval, then L 2 (J;A) represents the Banach space consisting 
of all quadratically integrable functions / : J —> X (in Bocher sense) with norm: ||/(t)||L 2 (j;X) := 

(fj \\f (t)\\xdt) . We denote by iL 1 (J; X) the space of all functions u £ L 2 (J;X) such that u', the 

weak derivative of u with respect to time variable, exists and belongs to L 2 (J; X), endowed with the 
/ \ 1/2 

norm ||u|| ff i( J;A ) = ( II u IIl2(j ; .y) + \\ u '\\l 2 (j-,x )) • For 1 < hj < n, we write Diu = du/dxi and 

Di.ju = d 2 u/dxidxj. For u £ H 1 (U) and f £ H 1 ( J; AT 1 (C/)), we also set the semi-norms \u\h^{u) '■= 
||Vu|| L 2 (c/) and \/\l 2 (j-h 1 (u)) '■= (fj For ease of notation, we do not always distinguish 

between the notation of u, u(t), u(t,x ) and Sometimes, the notation is not changed when a 

function defined on f 1 restricted to a subset. For the sake of brevity, we systematically use the expression 
A < B to indicate that A < CB for constant C that is independent of A and B. In some special cases, 
we may specify the used constants. 

2 Preliminary 

Let fl be a bounded domain with Lipschitz boundary. Let a and e denote the conductivity and 
permittivity distributions inside f2. We assume that a and e belong to L°°(f2). Biological tissues induce 
capacitive effects due to their cell membrane structures m- When they are exposed to electric pulses, 
the voltage potential it is a solution to the following time-dependent equation mm 

( —V-(a(x)Vu(t,x)+e(x)Vu'(t,x)) = f(t,x), (t, x) £ (0, T) x fl, 

^ u = 0, (t, x) £ (0, T) x dCl, (2.1) 

[ it(0, x) = uq , x £ ST, 

where uq is the initial voltage and T is the final observation time and / £ L 2 (]0, T\; 7L _1 (fl)) is the 
electric pulse. 

The goal of this work is to establish the well-posedness of the model system (12.11) and derive a fully 
discrete finite element scheme for the numerical solution of the system. Of our special interest is the 

case when the physical coefficients are discontinuous in fl, namely they may have large jumps across the 

interface between two different media, which is a common feature in applications, and the conductivity 
distribution <r(: r) does not need to be bounded below strictly positively. As far as we know, this is the 
first mathematical and numerical work on pulsed electrical fields in capacitive media. The main difficulty 
comes from the fact that (ED does not belong to the well-studied classes of time-dependent equations. 
Our results in this paper have potential applications in cell electrofusion and electroporation using eletric 
pulses [19] and in electrosensing (D- 

In this section, we first introduce some notions and preliminary results. For the sake of brevity, we 
write I =]0,T[, H = L 2 (fl), V = Bg(^) with its dual space V' = H~ 1 (Q) and X = H^fl) D iJ 2 (fl). 
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Clearly, V C H C V' is a triple of spaces (cf. [253 Chapter 1]), i.e., 

(1) the embeddings V C H C V' are dense and continuous; 

(2) {V ', V} forms an adjoint pair with duality product (■r)v'xv; 

(3) the duality product (-,-)v'xV satisfies 


{u,v)v'xv = (u,v), V u £ H, v £ V. 


We also introduce two bilinear forms ai(u,v) and a 2 (u, v) on V as follows: 


ai(u,v) = / a(x)X7u(x) ■ S7v(x)dx, 0 , 2 ( 11 , v) = / e(x)’Vu(x) ■ S7v(x)dx, u,v£V. (2-2) 

J n Jo, 

We first define the weak and strong solutions to the equation (12. 1 1) . We adapt the widely used notions 
of weak and strong solutions of parabolic equations (see, for instance, [2D]). 

Definition 2.1. Let uo £ V and f £ L 2 (I;V'). A function u £ lf 1 (/;y) is called a week solution of 
m if it(0) = u 0 and it satisfies the following weak formulation: 

ai(u(t,-),v) +a 2 (u'(t,-),v) = (f(t,-),v) v > xV (2.3) 

for all v £ Hq(Q) and a.e. t £ I. 

Definition 2.2. Let f £ L 2 (I ; H) and uo £ X. Then, a function u £ H l (I\ X) is called a strong solution 
of \2. 7]) if u( 0) = uq and the relation 

— V • ( o(x)S7u(t, x) + e(x)\7u'(t, x)) = f(t, x ) (2.4) 

holds for a.e. t £ I and a.e. x £ Cl. 

Remark 2.3. Let X be a Banach space. From [201 Proposition 7.1] we know that U 1 (/;X) (e C(I~X) 
continuously and 

sup||u(t)|| < ||u||j?i ( r ; x). (2-5) 

tei 

In particular, we have that u £ C(I ; V) for u £ R 1 (/; V). 

To prove the existence below, we will use the so-called “method of continuity”, whose key tool is the 
following lemma (e-g., 0). 

Lemma 2.4. Let X be a Banach space, Y a normed linear space, and Lq, L 1 two bounded linear operators 
from X to Y. For each X £ [0,1], set 

L\ = (1 — X)Lq + ALi, 


and suppose that there exists a constant C such that 


\\x\\x <C\\L x x\\ y , Vxel, A €[0,1]. 


Then L\ maps X onto Y if and only if Lq maps X onto Y. 


Let u be a function in a domain U C R ra , W <s U and the unit coordinate vector in the Xk direction. 
We define the difference quotient of u in the direction e& by 


D%u(x) 


u{x + hek) — u(x) 
h 


( 2 . 6 ) 


for x £ W and h £ R with 0 < \h\ < dist(W) dU). We will use the following lemma in the proof of 
Theorem 13.51 concerning the difference quotient of functions in Sobolev spaces (cf. [5] Lemma 7.23]). 
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Lemma 2.5. Suppose that u £ H l (U). Then for each W (e U, 

\\Dku\\ L 2 (w) <\\D k u\\ L 2 {u) , V h : 0 < \h\ < idist(W, dU). 

We end up with an analogue of [5] Lemma 7.24], and provide an outline of its proof. 

Lemma 2.6. Let u £ L 2 (/; L 2 (f/)), W <e U and suppose that there exists a positive constant K such 
that \\D%u\\ L 2 (!. L 2 ( W ^ < K for all 0 < \h\ < \dist{W,dU). Then \\D k u\\ L 2 ^. L 2 ^ w ^ < K. 

Proof. Banach-Alaoglu theorem implies that there exists a sequence {h m }^ =1 with h m —> 0 and a 
function v £ L 2 (/; L 2 (W)) such that \\v\\l 2 (i- L 2 (w)) < K, and for any ip £ C£°(W) and a £ Cg°(I), 


a(t)ipD k m u[t)dxdt —> 


i Jw 


a(t)ipv(t)dxdt 


i Jw 


as to —> oo . 


On the other hand, we have 


'k 

U Jw 

as to —> oo. Hence, we have 


a(t)pD^ m u{t)dxdt = — / / a(t)u(t)D k hm ipdxdt —> — / / a(t)u(t)Dk<pdxdt , 


i Jw 


i Jw 


/ / aft) (u{t)Dkip + v(t)p) dxdt. = 0. 

Ji Jw 


Using the arbitrariness of a and tp, we know for a.e. t £ I, v(t) = Dku(t) in weak sense, hence v = DkU 
in L 2 (I-L 2 (W)). □ 

Lemma 2.7. Let U C M™ be a domain and 1 < i < n. If u £ i7 1 (/; L 2 (U)), Diu' £ L 2 (I\ L 2 {U)) and 
Diu( 0) £ L 2 (U), then Diu £ L 2 (/;L 2 (C/)) and 

\\DM\l 2 {I-L 2 (U)) fz, \\^ > i u '\\L 2 (I-L 2 (U)) + ||Aw(0)||x,2((/). 

Proof. Since u £ L 2 (U)), we have 

u(t) = u(0) + f u'(s)ds, Vt£l. 

Jo 

By Fubini’s theorem, we know that for any <j> £ Cq°(U), 

[ u{t)Di(j)dx = f u(0)Di(f>dx + f f u\s)Di(j)dsdx = — f ( Diu( 0 ) + I DiU 1 {s)ds \ (fdx, 

Ju Ju Ju Jo Ju \ Jo J 

which implies that 

Di,u(t) = Diu{ 0) + / Diu'{s)ds. 

Jo 

This completes the proof. □ 


3 Existence and regularity 

We now introduce a basic assumption for the existence and uniqueness of weak solutions to (ED. 

(HI) a, e £ L°°(f2), and there exist two positive constants m and M such that 0 < <r{x) < M and 
to. < e(x) < M for a.e. x £ LI. 
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Let us recall that there exist two operators Ai,A2 : V —> V' associated with the bilinear forms ai(-,-) 
and a 2 (-,-), respectively, i.e., 

{AiU,v)v xV = ai(u,v), (A 2 u,v) v > xV = a 2 (u,v), u,v£V. 

>From [25] Theorem 1.24] we know that Ai is a bounded operator and satisfies the following estimate 

\\A lU \\ v , <M\\u\\ Vl V«eh, (3.1) 

and A2 is actually an isomorphism from V to V' and satisfies 

m||u||v < H^ullv' < M||u||v) Vueh. (3.2) 


3.1 Existence and uniqueness of weak solutions 

In this subsection, we prove the existence and uniqueness of the weak solutions to (12.11) . The first 
auxiliary result is the following a prior estimate, which lays the foundation for our subsequent existence 
and regularity results of weak solutions to m- 

Theorem 3.1. Let f £ L 2 (I;V'), uq £ V and u be the weak solution to (12.11) . Under the assumption 
(Al), we have 


/ 0 Jn 


\S7u’\ 2 dxdt + sup \\u{t)\\ 2 v < ||/||| 2 (J;V , 0 + |Kllt>, 


tei 


and 


IMIffhqv) ^ ll/IU 2 (/;W) + IMIw 

Proof. Choosing v = u' in (12.31) and integrating over (0,T), we obtain 


(3.3) 

(3.4) 


Jo J n 

>From this and the identity that 

r T r 


(a\7u(t) ■ Vu'(t) + e|Vi/(t)| 2 ) dxdt < / \\f(t)\\v'\\u'(t)\\vdt. 


cr V u(t) ■ Vu\t)dxdt = || v / o r Vu(T)||^ - || v / ^Vu(0)|| 


Hi 


Jo J n 


it follows that 


e\Vu'(t)\ 2 dxdt < / \\f(t)\\ V '\\u'(t)\\vdt +M\\u 0 \\ v . 


Jo Jn 

Using Young’s inequality, we have 


(3.5) 


(3.6) 


\\ u '\\l 2 (I-V) ^ ll/IU 2 (/ ; y') + I|mo||v- 

>From Lemma T2.71 and Remark 3.2, the desired results follow immediately. □ 

With estimate (13.41) in hand, we can prove the first existence result of (12.11) . 

Theorem 3.2. Let f £ L 2 (I;V') and uq £ V. Under the assumption (Al), equation (12.11) admits a 
unique weak solution. 

Proof. We first establish the results for uq = 0. The uniqueness is nothing but a direct consequence 
of Theorem o We use Lemma m to prove the existence. First, we construct a linear operator 
C : Hq(P, V) —> L 2 (I-,V') by setting 

(. Cu)(t) := Aiu(t) + A 2 u'(t), Vu £ Hq(I; V), 
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where ILq (/; V ) is defined by 


H 1 0 (I;V) = {ueH 1 (I-V);u( 0) = 0}. 

It is a closed subspace of the Banach space H 1 (I;V), since 7J 1 (/;y) (e C{I\V) continuously. >From 
m and (13.21) it follows that 

11 ^ 11 ^ 2(7 ^/) < M\\u\\ H i(iy), 

which implies that £ is well-defined and continuous. 

For each A G [0,1], we introduce a linear operator £\ : (/; V ) Lr{I ; V') as follows: 

£\u := A Cu + (1 — A)£ou, Vw G Hq{I] V), 

where we set £qu = —Au — Au' for u G iLg (/; V). Here (—A) is seen as an operator from V to V' (cf. [251 , 
Theorem 2.2]). More precisely, it is the operator associated with the bilinear form a(-,-) : V x V —> R, 
defined by 

a(u, v) := / \7u-S7vdx, Vu,uGH, 
in 

in a way such that ((—A)u,v)v'xv = a(u,v) for all u,v G V. In addition, (—A) : V —> V’ is an 
isomorphism. 

Let a\ = Act + (1 — A)xn and £\ = \e + (1 — A)xn for A G [0, 1]. Then the functions cr\ and £\ satisfy 
m! := min{m, 1} < E\(x) < M 1 := max{M, 1} for a.e. x G H, 


and 

0 < a\(x) < M' for a.e. iGf l. 

Then, for / G L 2 (I; V'), if C\u = / for some u G (/; V'), then u(0) = 0 and for a.e. t G I, u satisfies 
the following weak formulation: 

I (a\\7u(t) ■ Vu + £\\/u'(t) ■ X7v) dx = v)vxv, Vu G V. (3.7) 

Jn 

Thus, an application of Theorem 13.11 yields that there exists a positive constant C, depending only on 
to', M’ and T, such that 

ll u l|ff 1 (7;V r ) < C\\L\u\\ L 2(j v ,y 

In view of Lemma [2.41 it remains to prove that the mapping Cq : Hq(I; V) —> L 2 (I ; V') is onto. To this 
end, for an arbitrary / G L 2 (/; V'), we need to construct a function w G Hq(I; V ) such that for a.e. i G I, 

a(w(t),v)+a{w'(t),v) = {f(t),v) v >xv, Vv G V. (3.8) 

Let w(t) = J* e~ t+s h(s)ds for 1 G I, where h(s) = (—A) _1 /(s) for s G I. Since (—A) -1 : V' —> V is 
bounded, we have that h G L 2 (I;V) and hence w G iL 1 (/;’F). Moreover, a direct computation yields 
w(0) = 0 and w'(t) + w(t) = h(t) for i G I, which ensures that w satisfies (13.81) . Therefore, we can 
conclude that the method of continuity applies and Theorem 13 . 21 holds for uq = 0. 

For wo 7^ 0, we let w G U 1 (/;V r ) such that w(0) = uq and write f* = A\w + A 2 U)'■ Clearly, 
f* G L 2 (I\ V'). Then the proof above for uo = 0 confirms the existence of a unique function v G H 1 (/; V) 
such that u(0) = 0 and 

Aiv + A 2 v' = f - f*. 

Therefore, the function u := w + v is the desired weak solution. □ 
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3.2 Regularity of the solutions to the interface problem 

In this subsection we consider the regularity of the weak solution for (12. II) . which is important not only 
for its theoretical interest but also for the subsequent numerical analysis. Of our prime concern in this 
paper is the case when the coefficients a(x) and e(x ) are discontinuous, and the conductivity distribution 
u(x) in ED is unnecessary to be bounded below strictly positively. This feature is common to biological 
applications. Due to the (possibly sharp) jumps of cr(x) and e(x) across the medium interface, the solution 
to (12.11) does not expect a desired global regularity like H 2 (Ll), but it is shown in this section that this 
H 2 -regularity is true locally in each medium region fi* for i = 1, 2 . And such local U 2 -regularity is 
proved sufficient in the next section for us to establish the desired optimal convergence order for finite 
element approximations. As we are not aware of proofs of such local regularities for time-dependent 
PDEs with large jumps in coefficients in literature, even for standard parabolic equations, we will present 
a rigorous proof here for the non-standard time-dependent PDE (12.11) . We start with the introduction of 
some standard assumptions. 

(A2) n consists of two C 2 -subdomains and D 2 with Qi <£ 0, f l 2 '■= 

(A3) Ei := ejoi and eq := ojn^ are continuously differentiable in f \ (i = 1 , 2 ). 

The interface problem (12.11) is often complemented with the following physical interface conditions: 

[u(t)l =0 on I x T, ^ duft) + e ^ u =0 on I xT, ( 3 . 9 ) 

ou ou 

where T := dLl x is the interface, and [u(t)} := iti| r - u 2 \v, := er ± du ^ + <r 2 d gl^ + 

Ei d g 1 l }p +e 2 = 57 ^ on r. Here Ui stand for the restrictions of u to fand d/dui denotes the outer normal 
derivative with respect to flj, i = 1, 2. To deal with the interface problem, we introduce a Banach space 

y = {ueV-,u i eH 2 (n i ), * = 1 , 2 } 

with the norm 

IMIy = IMIv + ||Mi||ff 2 (ni ) + ||n 2 1|jT 2 (n 2 )3 Vn e y. 

Definition 3.3. Let f £ L 2 {I\H) and uo £ jV- A function u £ H 1 (/;jV) is called a strong solution of 
(ED with the jump conditions EH if u( 0 ) = uq and the relation 

— V • (cr(x)Vn(t, x) + E(ajV«'(f, x)) = f(t, x) (3.10) 

holds for a.e. t £ I and a.e. x £ fh (* = 1,2). 

Before proving the existence of a strong solution to the interface problem, we first establish the 
following result. 

Lemma 3.4. Let u he the weak solution of (|2-3D . Assume that f £ L 2 (I; H), uo £ y, u £ 1/ 1 (/; y), dil\ 
and dfl 2 are Lipschitz continuous. Then u is a strong solution for ED and (1X91) . 

Proof. We obtain, upon integration by parts, that for a.e. t £ /, 

f (-X7-(crX7u + eX7u')v — fv)dx= ( {aS7u ■ Vu + eVn • Vu — fv) dx, Vu £ .ffg(fh), (3.11) 
J Eli J£li 

which implies that 

—V • (cr(x)Vu(t, x) + e(x)Vu'(t, x)) = f(t, x) 

holds for a.e. t £ I and a.e. x £ fli (i = 1,2). It remains to show that the weak solution also satisfies 
the jump conditions (13.91) . By integration by parts we have for a.e. t £ 7, 


0= / (—V • (aVu + eVu') v — fv) dx 

J Uf22 

f f du du' 

= / (aVu ■ Vv + e\7u ■ \7v — fv) dx — / \a— —(-e—— \vdx, 
Jn Jr du du 


Vu £ V. 
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>From this and the definition of weak solutions it follows that 


/' du du' 

[a——\-£——]vdx = 0 1 \/v£V. 

J p uis ov 

The arbitrariness of v shows that u satisfies the second jump condition in (13.91) . The first condition in 
(13.91) is a direct consequence of the fact that u £ H 1 (I\V). This completes the proof. □ 

>From the lemma above, we know that the key point is to get the regularity of the weak solutions, 
which is the main subject of the following theorem. 

Theorem 3.5. Let f £ L 2 (I;H) and uq £ JL Under the assumptions (Hi), (A2) and (H3), the interface 
problem EH) and El admits a unique strong solution u, which satisfies 

~ ll/IU 2 (/ ; ff) + IKIly- 

Proof. >From Corollary 13.21 there exists a weak solution u £ H 1 (I; V ) to (12.11) . In view of Lemma EH 
and Theorem EH it suffices to show that u £ H 1 (I-,y) and 

\\ u \\ H^(I-y) ^ \\u\\m (I-V) + \\f\\ L 2 (I;H) + llwolly- 

The proof is divided into two parts. We only show that win! £ Ff 2 (fli), since the result that u\n 2 G 
H 2 (UL 2 ) can be proven in the same way. Henceforth we denote by C a generic constant that depends 
only on the cut-off functions, the final observation time T and the coefficients e and cr, and is always 
independent of the size of the difference parameter h in (12.61) . 

We first establish the interior regularity of the solution and its desired estimate. Let U <s Hi and 
choose a domain W such that U <s W <s Hi. We then select a cut-off function rj £ Cq°(W ) such that 
r/ = 1 on U and vanishes outside of W. Now let \h\ > 0 be small, and be the unit coordinate vector 
in Xk direction for k £ {1, • • • ,n}, and define a function v = —Df^lr^D^u!) (see (12.61) for the definition 
of -D£). Clearly, we know v(£) £ ffg(fli), hence also to V for t £ I. Now, letting a h (x) = a{x + hek), 
e h [x ) = e(x + hek) for x £ W, substituting this v into the left-hand side of (12.31) . and integrating it over 
/, we find that 


A := J (ai(u(t), v(t)) + a 2 (u'(t ), v(t))) dt 

= [ ! {Dl(aWu{t))-\l(p 2 Dlu\t))+D h k {£S/u'{t))^( V 2 D h k u'{t)))dxdt 

J i J n 

= [ f (E k t| 2 Vi) fe k t ( '(t) • Vfl^'(f) + <r' 1 r) 2 V^u'(f) ■ Vfi^(f)) (tedt 
JI Jw 

+ I I (2i 1 Dl£Dlu\t)\7u'{t)-S7r 1 + rj 2 Dle\7u'{i)-S7Dlu'{t) + 2£ h r 1 Dlu'{t)S7D } ]:u'{t)-\7i 1 )dxdt 

Ji Jw 

[ [ (2r)D k <jD k u' (f)V«(t) • V? 7 + rj 2 D k (jX7u(t) ■ X7D k u'(t) + 2a h gD k u' (t)V D k u(t) ■ Vry) dxdt 
Ji Jw 


n Jw 

=: (J)i + (*1)2 + (J )3 ■ 


We now estimate (J)i, (J)2, and (J)3 one by one. It is easy to see that 

= a h r 1 2 \S7D h k u{T)\ 2 dx- \ [ a h V 2 \\7D h k u(0)\ 2 dx + f f e h r, 2 \X7D h k u\t)\ 2 dxdt. (3.12) 
z Jw z Jw j I Jw 

We note that there exists a constant K > 0 such that \D k a(x)\ < K and \D k s(x)\ < K for all x £ W 
and 0 < \h\ < idist(W, 9f2!). Using Young’s inequality and Lemma 12.51 we obtain 


|(J)2 |<^/ f v 2 \VD^u'{t)\ 2 dxdt + C f [ \\7u'(t)\ 2 dxdt. 
5 J 1 Jw Ji J n 


(3.13) 


Similarly, we can derive 


|(J)s| < ^ f ( r/ 2 \VD^u'(t)\ 2 dxdt + 6 f ( r/ 2 \V D^u(t)\ 2 dxdt 

5 Jijw JiJw (3.14) 

+C (|Vu(f )| 2 + |Vu'(t)| 2 )dxdt, 

Ji Jn 

where <5 is a positive constant to be specified later. An interplay of Lemmas 12 . 51 and 12.71 implies that 

f y\\7D%u(t)\ 2 dx < C'( [ \\7D k u(0)\ 2 dx + f I \\7D k u'(s)\ 2 dxds), Vi 
Jw Jn Ji Jn 


£ I, 


with some constant C' > 0, whence (13.1411 ensures that 


2 m 


|(J)s| < / / r] 2 \S7D k u'(t)\ 2 dxdt + C / / (\Vu(t)\ 2 + \Vu'{t)\ 2 ) dxdt, 


i Jw 


i Jn 


(3.15) 


if S is chosen small enough, say <5 = m/(5TC). 

On the other hand, using Young’s inequality and Lemma 12.51 again, we deduce 


B := f f f(t)v(t)dxdt 
Ji Jn 

< V 2 \VD%u'(t)\ 2 dxdt + C J (|/(f )| 2 + |Vzt'(i)| 2 ) dxdt + ||u 0 ||;y^ . 

Since A = B, we combine (13.1211 with (13.161) to get 

~T [ [ V 2 \VD%u'(t)\ 2 dxdt 
0 JI Jw 

< t[ f V 2 \VDku'(t)\ 2 dxdt + C([ [ (\Vu(t)\ 2 + \Vu'(t)\ 2 + \f(t)\ 2 )dxdt + \\uo\\y), 
o Ji Jw Ji Jn 

which implies 


E / \\D h k D iU 'm 2 LHu) dt< [ (iK(t)lly + \\u{i)\\ 2 y + \\f{t)\\ 2 H ) dt + lluoll 2 ; 

< 1—1 JI JI 


for all k = 1,2, • • • , n and sufficiently small \h\ ^ 0. By applying Lemmas 12.61 and 12 . 71 we come to 


IMIffi (I\H2(U)) ^ ll/IU 2 (/;ff) + Ikllffi (I-V) + llwolly- (3.17) 

Next, we establish the boundary regularity and the desired estimate. We first use the standard 
argument to straighten out the boundary, i.e., flatting out the boundary by changing the coordinates 
near a boundary point (cf. [8] Chap. 6.2]). Given Xq £ SHi, there exists a ball B = B r (x o) with radius 
r and a C 2 -diffeormorphism 4/ : B —> ^l(B) C R" such that det|V'L| = 1, U' := ^(B) is an open set, 
d'(iJnfli) C R” and ^(.BflSfii) C 0R", where R" is the half-space in the new coordinates. Henceforth 
we write y = ^(rc) = (4'i(a;), • • • , 4 , n (*)) for x £ B. Then we have {y n > 0; y £ U’} = 4/(5 D Hi). Let 
$ = 4/ _1 , B + = B^{x o) D Hi, G = 4 '(i?r (xq)) and G + = ^>(B + ), then we can see G (s V and G + C G. 
We shall write DiW = dw/dyi for i = 1, • • ■ , n, and w(y) = u($(y)), f(y) = /($(?/)) for y £ U'. Now 
using the transformation function \&, the original equation on I x B can be transformed into an equation 
of the same form on / x U', i.e., for a.e. t £ /, 



&ijDiw(t)DjV + E £ijDi w '(t)Dj v )dy 

i,j=1 


[ f(t)vdy, 
Ju‘ 


Vv £ H]{U '), 


(3.18) 
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where the coefficients (j l:j ( y ) and ( y) are given by 


<Mj/) : = 4/(y) ~J2 £ (^(y))^r(^(y))^f-(^(y)) ( 3 - 19 ) 

r —1 r —1 

for 1 < i, j < n and y £ U'. It is not difficult to see that 

n n 

y £ ij(y)^j > H£l 2 > y > o, v (?/>£) g u' x M n , 

i,j=l j,i=l 


Choosing a domain W 7 such that G <s W 7 <e 17', we then select a cut-off function, which is still 
denoted by y, such that y = 1 on G and vanishes outside W'. Now let \h\ > 0 be small, and Bk be the 
unit coordinate vector in the yk direction for k £ {1, • • • , n— 1}. In the sequel, D £ stands for the difference 
quotient in the direction Bk■ We observe that there exists a constant K' > 0 such that \D l l '.d l ,j(y)\ < K' 
and \D%£ij(y)\ < K' for a.e. y £ W', all 0 < \h\ < |dist (W',dU') and 1 < i. j < n. Then, a natural 
variant of the reasoning leading to (13.161) shows that 


< 



+ C\\w(Q)\\ H 2(jj I _ ijU i + ), 


J dydt 
dydt + C 


“h 1 {U') + ll^Wllff 1 ([/') + \\Ht)\\L 2 (U'))dt 


where ||iy(0)|| J j2 (I71ut/ /_ ) := ||w(0)|| ff 2 (t7i) + ||w( 0 )|| H 2 ([/ ^ ) with U' + = U' D R" and U'_ = U'\ U' + . We 
can derive from the resulting inequality that 

n p 

J2 / \\D h k D lW '{t)\\l 2{G+) dt 

i =1 

~ J (\\ w '(.t)\\ 2 H 1 {u l ) + 11^(0II+ ll/Wlll^t/')) dt + ||w( 0 )|| ff 2 ( C/ /_ u£/ ^) 

for k = 1,- • • ,n — 1 and all sufficiently small \h\ yf 0, where we have also used the fact rj = 1 on G + . 
Using Lemma [2.01 we have 

y, ||Ajw|| ff i(/ ;i 2 ( G +)) < ||/||l 2 (J ; L 2 ( (7 /)) + ||w||iJi(/ ; JJi(C/')) + ll u; (0)||i/2((7^ u ;y' h ), (3.20) 

l<i,j<2n 

where Dijw = DiDjW. >From (13.181) we obtain upon integration by parts that for a.e. t £ /, 



for any tp £ G“(G + ). Noting that dy, and are both continuously differentiable in G + and the estimate 
(13.201) . the right-hand side of the equation above is well-defined and so we find that for a.e. t, £ /, the 
weak derivative of d nn D n w(t) + £ nn D n w'(t) with respect to y n exists and it satisfies 

- D n (a nn D n w(t) + £ nn D n w'{t)) = f(t) + y Di{£ijDjw'(t)) + y D^dijDjwit)). (3.21) 

l<i,j<2n l<i,j<2n 
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For the sake of brevity, we write g := i nn HD n w, where 

H(t,y) ■■= exp ( ann ^ t) (t,y) £ I x G + . 

£nn{y) 

It follows readily that H is strictly positive and ff S C J (/ x G + ). (13.211) ensures that for a.e. t £ /, 

-Dn(^P-)=f(t)+ Di&jDjw'it)) + Di(aijDj W (t)). (3.22) 

' l<i,j<2n l<i,j<2n 

A direct computation yields for a.e. t £ I, 

D n9 , (t) = ^W-g'(t)-H(t)(^m+ D z (iijD J w l (t))+ £ D^D^))) ■ (3-23) 

^ l<i,j<2n l<i,j<2n 

Since ||g , |U 2 (/;L 2 (G+)) "X, |M|Iir 1 (/;if 1 (G+))> we infer from (13.201) and (13.231) that 

I|Ai.9 , ||l 2 (/;L 2 (G+)) < ||/||l 2 (/;L 2 (!7')) + IM|,F/i(/ ; ,F/i((y')) + IM(0) \\h 2 (U'_UU^) ■ (3.24) 

As ||£> n ( 7 ( 0 )|| < |M( 0 )llff 2 (f 7 |)j an application of Lemma 12.71 yields 

\\D n g\\L 2 (i-L 2 (G+)) < II/I|l 2 (/;l 2 ((7')) + IMI m{i-,w{U')) + ||w?(0)|| h 2 (u'_ uuf) • (3.25) 

We can then conclude from (13.241) and (13.251) that 

\\Dn,nW\\ H 1 (I;L 2 (G+)) ^ II/II L 2 (I-L 2 (U')) + \\w\\H 1 (I-H 1 (U')) + 11(0) 11 .f? 2 (t/L UC/^_)- 
Combining this with estimate (13.201) . and transforming w back to u in the resulting inequality, we find 

IMI h i (i,h 2 (b+)) ^ ||/l|i 2 (/;ff) + IMI m{i-v) + IMIM (3.26) 

By choosing a finite set of balls {B ri / 2 {xi)}f =l such that it covers the boundary and then adding the 
estimates over these balls, we obtain the desired result. □ 

Using the standard arguments (cf. [TO] Theorem 3.2.1.2]) with some natural modifications and the 
estimates above, we can prove the following regularity result in a general convex domain. 

Theorem 3.6. Let f £ L 2 (I;H) and u o £ y. Assume that Q is a bounded and convex domain, fli <s U 
a C 2 subdomain, and that (Al) and (A3) hold. Then, the interface problem (12.11) and (13.91) admits a 
unique strong solution, which satisfies 

IMMc/fy) ~ \\f\\L 2 (i-H) + IKIM (3.27) 

3.3 Existence of a strong solution for smooth coefficients 

For the case with smooth coefficients, if we use 
(A4) dtt is C 2 and er, e £ C' 1 (U), 

instead of (A 2 ) and (A3), then we can obtain a better regularity result as follows, using the same argument 
as in the poof of Theorem 13.51 

Theorem 3.7. Let f £ L 2 (I;H) and uq £ X. Under assumptions (Al) and (Ad), the equation (12.11) 
admits a unique strong solution u, which satisfies the following estimate: 

IMI HUI-,X) \\f\\L 2 (I-H) + IMM- 

Remark 3.8. By the standard semigroup theory (cf. 1 25\). we can actually achieve a better estimate, 
i.e., under the assumptions of Theorem \3. 71 we have u £ C , 1 ([0, T]; X). That is, u is a classical solution. 


11 

















4 Finite element approximation and error estimates 

In this section we propose a fully discrete finite element scheme to approximate the solution of the 
interface problem (12.111 and (13.91) . and establish its optimal convergence under the minimum regularity 
assumptions on the given data. To do so, we first consider an auxiliary semi-discrete finite element scheme 
for the concerned interface problem and develop its optimal convergence, which will lead to the optimal 
convergence of the fully discrete scheme. 

Unless otherwise notified, we assume below that / £ L 2 (I; R) and u o £ y. For the sake of exposition, 
we further make the following assumptions: 

(A5) H is a convex polygon or polyhedron in R" with n = 2 or 3, and fix g SI is a domain with 
C 2 -boundary; 

(A6) The coefficients e and a are constants in each domain, namely, e = e, and cr = tji in flj, i = 1,2, 
where and eq are two positive constants. 

Clear, assumption (Al) is satisfied if (A6) holds. From Theorem 13.61 it follows that there exists a strong 
solution to the interface problem (12.11) and (13.91) . 

Remark 4.1. For the sake of exposition, we assume that Q is a convex polygon (if n = 2) or a convex 
polyhedral domain (if n = 3). The actual curved boundary can be treated in the same manner as we 
handle the interface T in our subsequent analysis of this section. 

We now introduce a triangulation of the domain Q. First we triangulate fli using a quasi-uniform 
mesh Tf with simplicial elements of size h, which form a polyhedral domain flq/j. The triangulation Tf 
is done such that all the boundary vertices of h lie on the boundary of Then we triangulate f l 2 
using a quasi-uniform mesh Tf with simplicial elements of size h, which form a polyhedral domain f 
The triangulation Tf is done such that all the vertices of the outer polyhedral boundary <9U are also the 
vertices of fl 2 ,h, while all the vertices on the inner boundary of Cl 2 ,h match the boundary vertices of f li t h- 
More precisely, the triangulation 7f satisfies the following conditions: 

(Tl) n = Ujc e7 - h K- 

(T2) if Ki,K 2 £ Th with K\ ^ If 2 , then either K\ D K 2 = 0 or K\ D if 2 is a common vertex, an edge or 
a face; 

(T3) for each if, all its vertex is completely contained in either Ui or U 2 . 

Now we define 14, to be the continuous piecewise linear finite element space on the triangulation 7h 
and V(( the closed subspace of 14 with its functions vanishing on the boundary dfl. Then, we study 
the approximation of piecewise smooth functions by finite elements in Vh- Clearly, the accuracy of 
this approximation depends on how well the mesh 7h resolve the interface T. Following the notation 
introduced in [15], we define, for A > 0 with A < min{dist(r, dfl), |}, a tubular neighborhood S\ of T by 

S\ := {x £ U; dist(a;,9r) < A}. 

Then, we decompose Th into three disjoint subsets Th = T^ U U %, where 

K = {K £ Tf, K c Qi\Sx}, * = 1,2, 

and T '■= Th\{Th U T((). Furthermore, we write Tf = {K gT*;K C Uj U S\}. Since T is of class C 2 , we 
know from mm that there exists A > 0 such that 

A = 0(h 2 ), (4.1) 

and T = Tf U Tff and T(f D T? = 0, provided that h is appropriately small. An important observation 
is that VL ih = U{if; If £ Tf U Tf}, i = 1,2, i.e., Tf = Tf U Tf. The notation S\ not only quantifies how 


12 


well the mesh Th resolves the interface, but it also allows us to use the lemma I4.3M.51 which were first 
established in m, in the subsequent analysis. 

We note that the evaluation of the entries of the stiffness matrix involving interface elements is not 
trivial in the three-dimensional case if the mesh is not aligned with the interface. So we shall adopt the 
following more convenient approximation bilinear forms cq^-, •) : V x V —> R: 



( 7 ,-Vu • \7vdx 


and 



£jVu ■ 'Vvdx. 


To approximate the problem in space optimally, we introduce the projection operator Qh : ynV —> V®. 
For each u G y, let f* = —EiAui in f i = 1,2, and g* = Clearly, f* £ H and g* £ L 2 (T). Then, 

we can define Qh : y D V —> V° by 


a 2 ,h(QhU,v h ) = ( f*,v h ) + ( g*,v h }, Vu G V°, 


where (■, •) denotes the scalar product in T 2 (r). We note that the right-hand side L(-) := (/*, •) + ( 5 *, •} 
is independent of h. Thus, we can follow the proof of ns Theorems 4.1 and 4.8], which mainly focuses 
on the case when g* = 0 , to obtain the following result. 

Lemma 4.2. We have 

a 2 (u,v h ) = a 2th (Q h u,v h ), Mv h G V°. (4.2) 

Moreover, for any u G y, the following error estimate holds: 


||m - Qhu\\ H + h\\u - Qhu\\v < h 2 \\u\\y. 


Now, we present some auxiliary results. For the difference between the bilinear form Gq(-,-) and its 
approximated bilinear form ■), we have the following result (cf. |T51 p. 27]). 

Lemma 4.3. Both au l (-,-) and a 2 ^(-,-) are bounded, and < 22 ,?i(t) coercive. Moreover, the form 
a^ h (u,v) := ai{u,v) — ai^h{u,v), i = 1,2, satisfies 

l a f ( M >' L ')I ^ Mi/ 1 (Sx)MiT 1 ('SA)- 

To estimate the energy-norm and the L 2 -norrn of a function over S\, we will frequently use the 
following result (cf. |15] Lemma 2.1 and Remark 4.2]). 

Lemma 4.4. For any u £V, we have 




(4.3) 

Moreover, for any u G y, 




U ^ A |u \\y 

(4.4) 

with | • |jyi(s x ) being the H 1 -semi norm. 




The following estimate is critical for proving our main result (cf. ]15] Lemma 4.5]). 
Lemma 4.5. There exists a positive constant p independent of h such that 




\j~-^\\ w h\\m(s^ h ), 


Vw h G V h . 
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4.1 Semi-discrete finite element approximation and error estimates 

We now consider an auxiliary semi-discrete finite element scheme for our concerned interface problem 
(ED and (13.911 and develop its optimal convergence, which will lead directly to the optimal convergence 
of the fully discrete scheme in section 4.2. 

Problem (P/,). Let Uh{ 0) = QhUo- Find Uh £ H 1 (P, V®) such that for a.e. t £ J, 

ai,h(uh(t),vh) + a 2 ,h(u' h (t),vh) = {f(t),vh)vxv, Vvh £ V°. (4.5) 


We first establish an auxiliary lemma, which will be used in the proof below. 
Lemma 4.6. Let f £ L 2 (P, V') and Uh be the solution to Problem (Ph)- We have 


\Wh\\m(I-V) $ ||/|| l 2 ( 7 ; 1 /') + \\QhUo\\v- 

Since the proof of this Lemma is the same as that of Theorem 13.11 we omit the details. With the 
results above, we are now in a position to prove the error estimate in the energy norm. 

Theorem 4»7« Lot u bo tho solutzoTz to tho zTztoT'fcico problem* m and ED and Uh the solution to 
Problem (Pk), then the following estimate holds: 

||« - u h \\m(i-v) < h (\\f\\Li(i-,H) + IKIly) • 


Proof. We first have the following decomposition: 


(IK*) - u h{t) \\v + IK(*) - u' h (t)\\v)dt 


< 


(IK*) - Qhu(t)\\y + || u'(t) - Q h u'(t)\\y)dt 


(4.6) 


+ |J (II Qhu{t) - u h {t)\\v + II Qhu'it) - u' h {t)\\y)dt 

=■■ (I)i + (I^- 


Using Lemma 14.21 and Theorem 13.61 we obtain 


(I)i ~ ^IMI H^i-y) ^ h (H/llff + IKIly) • 

It suffices to prove a similar estimate for (I) 2 . To this end, we first notice that the function w := Uh — QhU 
belongs to U 1 (/;V)( I ). In addition, using the identity that ( Qhu)'(t ) = Qhu'ff) for a.e. £ £ J and the 
definition of u and Uh, we find for a.e. t £ /, 

ai,h(w(t),v h ) + a 2 ,h{w'{t),v h ) = (F(t),Vh)v>xV, \/v h £ Vfi, 

where F(t) £ V' for t £ I, defined by 

( F(t),v)v x y := ai(u - QhU, v) + a 2 (u ' - Q h u',v) + a^{QhU,v) + a^iQhU 1 ,v), Vu £ V. 

Analogously to Lemma 14.61 we derive 


(fi2 = IMK(/;V) ~ II-F|Il 2 (/;V')- ( 4 - 7 ) 

Thus, it remains to estimate ||-F||l 2 (/;V')’ For t £ I and any v £ V, we use Lemma |4.31 to obtain 

\(F(t),v)vxv\ 

< (||u(£) - Q h u{f)\\v + \Q h u(t)\ H i {Sx ) + II u'(t) - Q h u\t)\\ v + \Q h u'(t)\ H i {Sx) ) ||u||y, 
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which, together with the estimates 


\QhU(t)\ H i(s x ) < |w(*)|h 1 (S a ) + | u(t) - Qhu[t)\ H i(s x ), 


and 

\Qhu'{t)\ H i(Sx) ^ + I u '(t) ~ Qh,u'(t)\ H i( Sx ), 


implies that 


\\F\\l^{I-V') < II QhU - u\\ H i(I,V) + |M| ff i(/ ;ff i(s A )). 
Now Lemmas 14.21 and 14.41 together with Theorem 13.61 yield 


ll^lli 2 (/;y') ^ (h + VX) (\\f\\L 2 (i-H) + ll M o||y) ■ 
>From this, (HD and (14.71) . the desired estimate for (I) 2 is established. 

Now, we are in a position to prove the L 2 -estimate. 

Theorem 4.8. We have the following estimate in L 2 -norm: 


□ 


II u - Uh\\L^(I-H) < h 2 (\\f\\L 2 (I-H) + ||Mo||y) ■ 

Proof. For the duality argument, we define w £ H l (I\ V) and Wh £ H 1 ^', V such that for a.e. t £ I, 


ai(w(t), v) - a 2 (w'(t), v ) = ( u{t) - u h {t ), v), Vv £ V, 
ai(wh(t), v) - a 2 {w' h (t),v) = ( u(t) - u h (t),v), Vv £ V^°, 

which satisfies w(T) = Wh{T) = 0. That is, w*{t) := w(T — t ) is the weak solution of (12.11) with initial 
value u>*(0) = 0 and / replaced by u — Uh- Then Theorem 13.61 implies that 


IMI m(i-y) < ||u - u h \\ L 2 (i-h)- ( 4 - 8 ) 

Using the same argument employed in Theorem 14.71 with a natural modification, we find that 


||w — whWmii-v)) ^ h\\u — u)i\\l 2 (i-,h)- 


(4.9) 


By integration by parts with respect to the time variable, identity (14.21) and taking advantage of the 
Galerkin orthogonality for w — Wh and e := u — Uh, we know that 


and for a.e. t £ /, 


(a 1 (e,w h ) - a 2 (e, w' h ))dt 

(ai(e, w h ) + a 2 (e', w h ))dt + a 2 (w(0) - Q h { 0), w h {0)) 

= / (-af (u h , w h ) ~ a, 2 (u' h ,w h )) dt + af(<2/,u(0), w h (0)), 
Jo 

ai(w(t) - Wh(t),v) - a 2 {w'{t) - w' h (t),v) =0, Mv £ V°. 


(4.10) 


(4.11) 
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Applying (14.101) and (14.111) and integrating by parts with respect to time variable, we obtain 


INI 


L 2 (I-H) 


(ai(e, w) — a 2 (e, w '))dt 


/ (ai(e, w — w h ) - a 2 (e,w’ - w' h ) + ai{e,w h ) - a 2 (e,w' h ))dt 
Jo t 

/ (ai(u - Qh,u,w - w h ) - a 2 (u - Qhu,w' - w' h ))dt 
Jo ... 


- / (of {u h ,w h ) + af (u' h ,Wh))dt-\-a 2 {Qhu{0),w h (0)) 

Jo 

= / (ai(u - QhU,w - w h ) + a 2 (u' - Q h u',w - w h ))dt 

Jo T 

- j (ai(u h ,Wh) + a2(u' h ,Wh))dt 

+a 2 (Q h u(0), w h (0)) + a 2 {u{ 0 ) - Q h u( 0 ), w(0) - w h { 0 )) 

=: (II) 1 + (II) 2 + (II)3. 

By the Cauchy-Schwarz inequality it follows that 

(II) 1 < \\u - Q h u\\ L 2 {I . v) \\w - w h \\ L 2 (iy) + II u' - Qhu'\\ L 2 {I . v) \\w - w h \\ L 2 {I . v) 
Applying Lemma 14.21 the regularity estimate (14.81) and Theorem 14.71 we have 

(H) 1. < h2 (ll/IU 2 (/;ff) + ||«o|| y) ||e||z,2(j ;fl -). 

By Lemma 14.31 and the Cauchy-Schwarz inequality, 


(4.12) 


(H)2 ~ (l u h|.L 2 (J;.ff 1 (Sx)) + \ u 'h\L 2 (I-,H 1 (Sx)))\ w h\L 2 (I;H 1 (Sx))- 

Before we further estimate (II) 2 , we first bound uh, w and Wh in H 1 (S\)). Applying Lemma 14.41 

and Theorem ro we have 

||w/i||j?i(J;ffi(S x )) ^ ll e llff 1 (7;i7 1 (S'x)) + l|w||i/i(/ ; ffi(S A )) < (a/A + /i)IMIff H^y)- ( 4 - 13 ) 

On the other hand, using Lemma 14.41 and the regularity estimate (14.81) . it follows that 

Using Lemmas 14.21 and 14.51 the regularity estimate (14.81) . (14.91) and the condition 2A < h, we have 

|Nh||ffl(/;ffl(S A )) < V\h~l ||W/ l || ff l(/ ;ff l( S(iA )) 

< V\h~i (||w- WhWHi^H^S^)) + IN||ff 1( J;i7i(S MA ))) (4.14) 

< V\h~2 (j|uJ - W/ l ||iji(/;ifi(s MA )) + ^ 2 IMlL 2 (/ ; y)) 

% v / A||e|| L 2 ( /;ff ). 

Now, (14.131) . (14.141) and Theorem 13.61 yield 

( 11)2 ~ + Vxh) (||/||l 2 (7;J7) + INolly) ||e||L2(/;j7). (4-15) 

To bound (II) 3 , we first need to estimate |Q^rt(0)|^i(5 A ) and |u)^(0)|//i(5 A ). To this end, applying Lemmas 
14.21 and 14.41 we find that 

|QhW(0)|fl-i(5 A ) < \QhU(0) - w(0)|ff 1 (S A ) + l u (0)|ff 1 ('S , x) ~ (>/A + /l)||tt(0)||y. 
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On the other hand, from (12.51) it readily follows that 


11 Wfc(0)l I h 1 (s x ) < SU P ll w; /i( i )llff 1 (SA) < 
tei 

Similarly, we have ||w(0) — w/ l (0)||y < ||it? — Wh\\m(iy)- We thus find 

11^(0)IIh^Sa) ~ ~ v / A||e|| i 2 (/ . ff ). 


Summarizing the above estimates, we get 


(H) 3 ^ IQ^(0)|jji ( s a) |^(0)| h i (i s a) + ||<5hu(0) - u(0)||^||w(0) - tu h (0)||v (4.16) 

< (A + hVX + h 2 )||u(0)|| 3 ;||e|| i 2(7. ff ). 

Taking (14.11) . (14.121) . (I4.15D . and (14.161) into consideration, we can conclude the desired estimate. □ 


4.2 Fully discrete finite element scheme and error estimates 

In this subsection, we are going to formulate a fully discrete finite element scheme to approximate 
the solution to the interface problem (12.11) and (13.91) . We shall use the backward Euler scheme for the 
time discretization. Let us start with dividing the time interval I into N equally spaced subintervals and 
using the following nodal points: 

0 = t° <t 1 < • • • < t N = T, 

where t n = nr for n = 0,1, • • • ,N and r = T/N. For any given discrete time sequence {u rl }))' =0 in V 
and a function g(x, t ) which is continuous with respect to t, we can define 


d T w n = 


w — w 


71—1 


r = ~f 9 (;s)ds , n = 1,- ■ 


,N. 


Now, we propose a fully discrete finite element scheme to approximate the solution to the interface 
problem (12.11) and (13.91) . 

Problem (Ph,r)- Let = Qh^o- For each n = 1, 2, ■ • • , N, find uJJ G V f ® such that 


ai, h (v%,v h ) + a 2 ,h(d T v%,v h ) = (f n ,v h ), Vv h G V°. 


(4.17) 


For a discrete sequence defined in Problem (Ph.r), we can introduce a piecewise constant 

function in time by 

u h , T (;t)=u%(-), Vt G t n ], n= 1,2, ■■•,1V. (4.18) 

Then, we say that Uh,r is a solution of Problem (Ph.r), which is a fully discrete approximation of the 
solution to the interface problem (12.11) and (13.91) . In order to compute the error between Uh, T and u, it 
suffices to establish the error between Uh, T and Uh, i.e. the solution of the semi-discrete scheme (14.51) . 
since the error between Uh and u has been studied in Section 4.1. To this end, we need the following 
auxiliary result. 

Lemma 4.9. Let {Fn}0Li be a time discrete sequence lying in V' and u>° = 0. There exists a unique 
sequence such that for n = 1,2, • • ■ , N, 

ai, h (w%,Vh) + a 2 ,h(d r 'w^,v h ) = {F n ,v) v > xV , \/v h G Vfi . (4.19) 


Moreover, the sequence {uih}n=i das the following stability estimate: 


max 

l<n<N 


N 

\M<r^\\F n \f v , 

n—1 


(4.20) 
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Proof. The existence and uniqueness follows immediately from the Lax-Milgram theorem. Taking Vh = 
2rdw r f in (14.191) . then using the relation 


2Ta 1 , h {Wh,d T Wh) = ai th {Wh,Wh)-ax th {w'll X )+T 2 ai , h {d T wZ, d T w%), Vn= 1,-- 

and the coercivity of 0 , 2 ,h (■, •), we obtain that 

2mT\\d T w%\\lr + ai , h {wZ,w%) - ai,/i(w^ -1 , w^ -1 ) < 2T\\F n \\ v/ \\d T w%\\ v , Vn = 1,2,--- 

Adding the above inequalities from n = 1 to n = TV, and using the Cauchy inequality, one has 

JV N 


■ ,iV, 

,1V. 


n =1 n=l 

Now the desired estimate follows immediately from the inequality 

N 

IKiry<Tr^R<|| 2 y , V 1 < n < TV. 

n= 1 


□ 


>From the lemma above we notice that Problem (Ph.r) always admits a unique solution. 

Lemma 4.10. Let Uh,r and Uh be the solution of Problem (Ph,r) and Problem (Ph), respectively. Under 
the assumption that f £ H l (I]H), the following estimate holds: 

IK ~ Ur,h\\L 2 (I-V) ^ T (||/ , ||l 2 (I;J?) + l|/l|i 2 (/;i/) + KHy) , 

Proof. We first define a piecewise constant function in time such that u £ T (0) = QhUo and 

u* htT (;t)=v %(•), Vi G n = l,2,---,lV. 

Using Lemmas 14.21 and 14.61 it follows readily that 

K - U* h ^\\ L 2 (I- V ) < T\\u h \\ H i (J ; y) < ^~( II /II Z/ 2 (/;//) + Klly)- (4-21) 

Integrating (TQ>1) over (f" 1 ,t n ) and dividing both sides by r, we have for n = 1,2, • • • ,7V, 

ai, h (uh,v h )+ a 2 ,h.{d T Uh,Vh) = {f n ,v h ), Vv h £ V°. (4.22) 

Subtracting both sides of (14.221) above from those of (14.171) . we can rewrite the resulting equation as 

ai,h( u h ~ Uh’Vh) + a 2 ,h(d T {ul - Uh),v h ) = (/" - ~f ,v h ) + a ljh (u% - u h ,v h ), Vv h £ V°. 

The right-hand side of the equation above defines a functional on V for each n = 1, 2, • • ■ ,N. Indeed, we 
have for n = 1,2, ■ ■ • , IV, 


|(/ n ~T,v) + a lth (u% - u h ,v)\ < (II/" - 7"||ff + IK -u h \\v) IH|y, V v £ V 


by using Poincare’s inequality. Therefore we can apply Lemma 14.91 to obtain 

N 

\\ u h,T ~~ u hA\ 2 L 2 (I-V) = E T \\ u h ~ u h\\v ^ ^ IK ~ u h\\v 


n—1 
N 

< r 

n—l 


E 


-f n \\H + \\u n h -U h \\ 2 V 


~ (II /^ IIL 2 (/;//) + \\ u h\\ 2 m(I-V)) 

< T2 (ll/ / Hl ( / ; ff ) + ll/lli ( J;J? ) + IKIly)- 
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Now the desired result follows from the previous estimate, (|4.21l) and the following triangular inequality 


\\ u h ~ Uh, T \\L 2 (I-,V) < II u h - «£, T IU 2 (/;V0 + || U* h T - Uh,T ||L 2 (/; V) ■ 


□ 


>From Lemma T4. 101 and Theorems 14.71 and 14.81 the following result follows immediately. 

Theorem 4.11. Let u be the solution to the interface problem (12.11) and (13.91) and Uf l T the solution to 

Problem (P h,r). Under the assumption of Lemma \f.lO\ the following estimates hold: 

||m - || L 2 (I;V) ^ (T + h ) (||/||_L2( /;ff ) + \\f'\\ L 2(i-,H) + Ikolly), 

||m - u h , T \\ L 2(i. H ) < (t + h 2 ) (||/||l2( /;H) + ||/'||l 2 (/;H ) + ||M 0 ||y)- 
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